A theory of 7r/2 superconducting Josephson junctions 
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Abstract 



1-^ ' We consider theoretically a Josephson junction with a superconducting criti- 

a' 

O ' cal current density which has a random sign along the junction's surface. We 



show that the ground state of the junction corresponds to the phase difference 



^ ' equal to ^. Such a situation can take place in superconductor- ferromagnet- 

Q I superconductor junctions. 
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The superfluid density connects supercurrent with superfluid velocity in superconductors. 
In the framework of BCS theory it is a positive quantity. Thus the ground state of a 
superconductor corresponds to zero superfluid velocity (See for example t^'). The question 
whether and under what circumstances the superfluid density can be negative has been 
discussed in many theoretical papers, a). It has been shown t^^^l that if the distribution 
function of quasiparticles in superconductors is non-equilibrium the superfluid density can be 
negative. Negative Josephson coupling has been observed experimentally in nonequilibrium 
superconductor-normal metal-superconductor junction ^^'. b). In the presence of magnetic 
impurities or resonance states with strong electron-electron correlation inside the insulator a 
superconductor-insulator-superconductor junction could have negative Josephson coupling 
[7-10]^ c). In case of disordered superconductor-ferromagnet-superconductor junctions the 
Josephson coupling averaged over samples has oscillating sign and decays exponentially with 
the ferromagnetic layer thickness L t^^l . Mesoscopic fluctuations of the coupling, however, 
do not decay exponentially with L. Therefore at large L the junctions have couplings with 
random sign l^^l. d). At last, a d-superconductor-normal metal-d-superconductor junction 
can have arbitrary sign depending on the relative orientation of the (superconducting) order 
parameters in the superconductors I^^'^^l. The fact that the superfluid density is negative 
indicates that the state with zero phase gradient is unstable. The ground state of a Josephson 
junction with negative critical current corresponds to a phase difference on the junction 
which equals vr. 

In real situation, however, the sign of the superfluid density can fluctuate from point to 
point inside the sample. In this paper we consider a Josephson junction, with a coupling 
J{p) fluctuating from point to point along the junction's surface and with a random sign. 
Here p is a coordinate along the junction's surface. We show that in this case the ground 
state of the system corresponds to a phase difference of the order parameter of ^. 

The total energy of the system can be written as 

E = Ej + '^ fd'r(vs{r)y (1) 



Here m is the electron mass, Vs{r) = — V0(r) and 0(r) are the superconducting velocity 
and the phase of the superconducting order parameter in the bulk of the superconductors, 
{h = 1), Hg is the density of superconducting electrons, which is assumed to be the same in 
both superconductors, Ej is the Josephson energy of the junction which can be written as 

Ej = - [cPpJ{p)cos{x{p)\ (2) 



x(p) = (f){z = +0, p) — (f){z = —0, p) is the jump of the phase of superconducting order 
parameter at the junction's surface z = 0. Here z is the coordinate perpendicular to the 
junction's surface. In Eqs.1,2 we neglected the thickness of the junction. The second term 
in Eq.l corresponds to the energy associated with supercurrents in the bulk of the supercon- 
ductors. We neglect the magnetic field energy in the bulk of the superconductors assuming 
that J{p) is small and the Josephson penetration length of the magnetic field along the 
junction's boundary is large. 

In the case, when J(p) has a random sign, the ground state of the system corresponds 
to nonuniform distributions of Vs{r) ^ and x{p) 7^ 0- The current conservation law in the 
bulk of superconductors and at the junction's boundary gives 

A0(r) = (3) 

d(j){p,z = +0) d(f){p,z = -0) ^ jt ^ ■ t t ^^ u\ 
dz = dz = n/^P^ ""(^(^)) (') 

respectively. Solutions of Eqs.3,4 determine the extrema of the total energy Eqs.l and 2. 

At given (x(p)) the minimum of Eq.l corresponds to some distribution of 5x{p)-, which is a 

random sample specific quantity. If J{p) is small enough, then (5x(p) = x(p) ~ (x(p)) ^ 1- 

Here the bracket () stands for averaging over the junction's surface. 

Substituting x{p) ~ ixip)) ^^ l^ft hand side of Eq.4 and solving Eqs.3-4 we obtain an 

expression for the minimum of the total junction's energy per unite area, 

Eeff = -{J{p)) cos(x) - Jeff sin^ix) (5) 

as a function of (x(p)), where 



■Jeff = -w- / dpidp2^— 6 

Here S is the area of the junction. For < 2 Jg// < (J) the junction's energy Eeff has a 
minimum for (x) = 0, while at 2Je// > (J), the minimum corresponds to a nonzero phase 
difference 

(x)'"^" = arccos(-^), (7) 

which in the case of Jeff ^ (J) gives (x)™*" = f • The derivation presented above is 
similar to the case of a ferromagnet-normal metal-ferromagnet junction with a random in 
sign exchange energy between the ferromagnets ^^^' 

Below we consider in more details the case of a superconductor- disordered ferromag- 
netic metal-superconductor junction shown in Fig.l. We will describe the system by the 
Hamiltonian 

H = Hbcs + Ht + Hp- (8) 

where Hbcs is the BCS Hamiltonian of superconducting leads, 

HT = t J2 I lldr^('^t{r^,a)^F{r,a) + h.c). (9) 

is the Hamiltonian describing the tunneling between the two superconductors labeled by 
indexes 1,2, and the ferromagnet labeled by the index F\ and a; is a spin index. Furthemore, 

Hf = Ho + Ha (10) 

is the Hamiltonian of the ferromagnet. The integration in Eq.9 is taken over the 
superconductor-ferromagnet surfaces. The effective magnetic field H models the ferromag- 
netic exchange spin splitting, cr are Pauli matrices; Hq is the hamiltonian of free electron gas 
in a random potential U{r). We assume the following correlation properties of the random 
potential : < U{r) >= and < U{r)U{r') >= ^wi^t ■ Here r is the electron elastic mean 
free time, uq is electron density of states at the Fermi level. 



Let us introduce a nonlocal Josephson coupling density J(p, p') as 

Ej= j dpdp'Jip, p') cos(0(p) - (Pip')) (11) 

where p and p' are coordinates on the surfaces between the ferromagnet and the two super- 
conductors. The integration in Eq.ll is taken over the surfaces S and S' of the supercon- 
ductors. 

To lowest order in t we have 

Jip,p') = 

T f [ 1 

= ij^t'^ drdr' F+_„(e„; p; r)GaA^n] r, p')^^ ,-«(£«; p'] r')G-a~ai~en; r', p) + h.c. (12) 



en, a 



Here F^_^(e„; rir[) and FQ,^_Q.(e„; rir[) are anomalous Green functions of the superconduct- 
ing leads. Ga,a{^n', fi, ^2) are Green functions of electron in ferromagnets; e^ = (2n + l)7rT 
is the fermionic Matsubara frequency. A derivation identical to Eq.6 gives 



Jeff = -^ / dpdpidp'dp\ 

O Tic J 



2tx\p — p' 
im ( g \V D \2 f 



TTUs VSttz/o-D/ \2-k'^L^J I 



Ht)- (13) 



Here L is the ferromagnet's thickness, ^ = JD/A is the superconducting coherence length 
and A is the modulus of the superconducting order parameter, g is the dimensionless conduc- 
tance of superconductor-ferromagnet boundary per unit area, and D is the electron diffusion 
constant, which is assumed to be the same in the superconductors and in the ferromagnet. 
Diagrams which contribute to Eq.l3 are shown in Fig.l. We use the conventional diagram 
technique for the averaging over scattering potential configurations '^^^ and retain the symbol 
() to indicate the averaging over configurations of the scattering potential. 

In the case of no electron reflection from the superconductor-ferromagnet boundary we 
can estimate J^ff by taking into account that the mesoscopic fluctuations of the critical 
current of a SNS junction of the area L^ is of order of -p- '^^1. Fluctuations of the Josephson 
coupling associated with different areas of the junction's surface separated by distances 
larger than L are uncorrelated. As a result, we estimate 



m / D^"^ 



It follows from Eq.l2 that Jeff decays with L only as a power law. On the other hand, 



(J(p)) decreases exponentially at L > J D/H t^^l and we have |(J(p))| ^ Jeff- 

In conclusion we would like to mention that the state discussed above of the 
superconductor- ferromagnet-superconductor junction is similar to Fulde-Ferrell-Larkin- 
Ovchinnikov state of bulk superconductors I^^'^^l, when in the presence of the spin split- 
ting energy the superconducting ground state corresponds to a nonuniform distribution of 
the order parameter. After averaging over the configurations of the impurity potential the 
order parameter in such a state is significantly suppressed by static disorder. Mesoscopic 
fiuctuations of the order parameter in such a state can however survive the disorder. 
Authors kindly acknowledge discussions with P. Fulde. 
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FIG. 1. A schematic picture of the superconductor (S)-ferromagnet (F)- superconductor (S) 
junction and the diagram for calculation of Eq.l3. Thin solid lines inside the ferromagnet cor- 
respond to electronic Green functions G{p,pi), dashed lines correspond to the elastic scattering, 
thin solid lines inside the superconductors correspond to the anomalous Green function F{p,p'). 



